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Abstract 

The relationship between pion and kaon space-time freeze-out distributions 
and the HBT radius parameters in high-energy nucleus-nucleus collisions is 
investigated. We show that the HBT radius parameters in general do not 



(^ [ reflect the R.M.S. deviations of the single particle production points. In- 

\-^ ' stead, the HBT radius parameters are most closely related to the curvature 

\^ ' of the two-particle space-time relative position distribution at the origin. We 

^^ ■ support our arguments by studies with a dynamical model (RQMD 2.4). 
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I. INTRODUCTION 

Two-particle intensity interferometry (HBT) has been used extensively in nucleus-nucleus 
collisions to measure the space-time extent of the particle emitting region [jl[ . Because of the 
complicated dynamics of high-energy nucleus-nucleus collisions, the HBT radii, defined as 
parameters of a Gaussian fit to the correlation function, do not measure the full size of the 
particle emitting region directly. Instead, the HBT radius parameters reflect an interplay 
between the geometric system size, the expansion of the system, the duration of particle 
emission, and the contribution of particles from resonance decays. It is hoped that the two- 
particle correlation functions can be parameterized in such a way that the various factors 
that influence the measured HBT radii can be understood individually. 

The two-particle correlation function is typically measured as a function of momentum 
difference q = pi — P2 and parameterized in term of radius parameters Ri and the so-called 
haoticity parameter A:[] 

C2(?) = l + Ae-''?^?. (1.1) 

For a static Gaussian pion source, the HBT radius parameters Ri directly correspond to the 
variance of the particle spatial distribution at freeze-out |l[] , 

R^ = RMSl^_y^,) = a\x, - y,t), (1.2) 

where Vi is the particle velocity and Ri are the HBT radius parameters obtained from a fit 
to the correlation function. For a non-Gaussian source, this simple relationship no longer 
holds. In this paper, we explore the meaning of the Gaussian radius parameters in the case 
of a non-Gaussian particle source and use a dynamical model to investigate the breakdown 



of equation 1.2 



II. BASIC RELATIONS 

Note that in general case of the particle source function the relation ( |1.2| ) cannot be true. 
A simple example would be the source with added a few particles emitted at very large 
distances. The addition of such particles can significantly change the freeze-out distribution 
R.M.S value and at the same time should not change the correlation function (sensitive only 
to pairs with close production points). Below we argue that the correlation function, and 
in particular the HBT radii derived from the Gaussian fit to the correlation function are 
indeed sensitive only to the distribution of close pairs. Hence, they can be quite different 
from the R.M.S. values of particle freeze-out distribution. 

Let us show that the Gaussian fit to the correlation function yields the value of {qf) = 
{2R^)~^ where the average is taken over the function C2{q) — 1. Consider a Gaussian fit to a 
one-dimensional correlation function using the maximum likelihood method. The normalized 
fit function is 



^ For simplicity, we omit in this parameterization the possible cross-terms [Q, which are not 
important for our study of particle production at mid-rapidity and in central collisions. 







The likelihood functional to be minimized is then 




InL = lnn/(?.) = InYli^Re-'^'''') = y:(ln4= + Ini?- 
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The likelihood equation becomes, 




ainL >r-/ 1 2] 




Solving this equation yields, 




1 _S^?'_/.2^ 





(2.1) 



(2.2) 



(2.3) 



2R^ N — ^'-'^ 

Hence, irrespective of the shape of the measured correlation function, the fitted Gaussian 
HBT radii measure the value of (g^). 

{q^) can also be related to the shape of the source function (time integrated relative 
distance distribution) S'(Ar), where Ar is the distance between particle production points 
Ti = Xi — Vit. Using this function, the correlation function in momentum space can be 
Fourier transformed into position space 0,H: 

C(g) - 1 = / e-^'i^'"5(Ar)dAr (2.5) 

5(Ar) = ^ / e'''''{C{q) - l)dq. (2.6) 

Differentiating the source function twice yields, 

1 d^S{/\Y) 



S{/\y) dArf 



W, (2.7) 



Ar=0 



where the average is taken over the correlation function, C{q) — 1. Eqns. |2.4| and |2.7| show 



that the HBT radius parameters are directly related to the second derivative of the source 
function of relative production points near Ar = 0. This is true irregardless of the shape of 
the correlation function. 

Note that the lambda parameter in the fit function (Eq. |1.1|) in a maximum likelihood 
approach is defined by the normalization 

/ ciqAe-^?^? = f{Ciq) - l)ciq. (2.8) 



The r.h.s. of this expression can be determined from equation |2.6| : 

J{C{q)-l)dci={2nfS{0). (2.9) 

It follows then that 

A = 2\^^^S{0)RiR2R3, (2.10) 

where S{0) is the (true) source function at the origin, and i?i, i?2, and R3 are the (fitted) 
radius parameters. 



III. MODEL CALCULATIONS 

To understand the relationship between the fitted HBT radii and the source function of 
relative production points and check how well the above equations work for a realistic particle 
source, the event generator RQMD (v2.4) [Q is employed. This model has been shown to 
reproduce experimental HBT data reasonably well 0. We focus on particle production 
in midrapidity region {\y\ < 0.5) in central (6 < 3fm) Au+Au collisions at RHIC energy 
(^5 = 200 GeV/nucleon). 

Figure |l| shows single particle freeze-out position distribution and the two-particle relative 
position distribution for pions. The single particle freeze-out position difference distribution 
has long non-Gaussian tails. In addition, the shape is not well described by a Gaussian 
at r = 0. The two-particle position difference distribution still shows large non-Gaussian 
tails. The distribution, however, is well described as a Gaussian near Ar = 0. Fitting this 
distribution to a Gaussian (also shown in the Figure) yields a Gaussian width a that is 
substantially smaller than the R.M.S. deviation of the distribution. 

The two-particle correlation functions are calculated as a function of pt for mid-rapidity 
{\y\ < 0.5) pions and kaons assuming plane wave propagation for the outgoing particles. 
The three-dimensional correlation function is then fit using the Bertsch-Pratt out-side-long 
parameterization, 

C2{qo, Qs, qi) = l + Ae-^"^"-"^'^^-"?^?, (3.1) 

where qi is the momentum difference along the beam axis, qo is the momentum difference 
perpendicular to the beam axis and parallel to the total transverse momentum of the pair, 
and qs is the momentum difference perpendicular to the beam axis and perpendicular to the 
total transverse momentum of the pair. 

Figure ^ shows the fitted HBT radius parameters as a function of transverse momentum. 
Also plotted are the R.M.S. deviations of the single particle freeze-out distributions and 
the fitted Gaussian widths of the two-particle position difference distributions. The R.M.S. 
deviations are substantially larger than the fitted HBT radii in all cases, while the fitted 
Gaussian widths of the two-particle relative position difference distributions are very close 
to the fitted HBT radii. The single-particle R.M.S. deviation differs from the fitted HBT 
radii most at low transverse momentum. This deviation arises from the long non-Gaussian 
tails in the freeze-out position distribution due to pions from resonance decay and particular 
expansion dynamics of the system. 

The width parameter (sigma) of the Gaussian fit to the two-particle difference distribu- 
tion near the origin (in our calculations we perform the fit in the region of ±0.5 R.M.S.) is a 
very good estimate of the second derivative of the relative source function at Ar = 0. Hence, 
the observed good agreement between HBT radii (fit parameters to the correlation function) 
and widths of the Gaussians fitted to the two-particle position difference distribution near 
the origin confirm relationship shown in eqns. |2.4| and |2.7| . Note a simple interpretation of 
the lambda parameter following from Fig. 1 (two-particle position difference distributions). 
It is just a product of the ratios of area under the Gaussian fit and the total number of pairs 
used for the distribution. 

Figure ^ is similar to figure || except that the particles used are kaons. In this case, 
the R.M.S. deviations of the single-particle freeze-out distributions and the fitted Gaussian 



width of the two-particle relative distance distribution are very similar in the transverse 
directions. Since fewer kaons come from resonance decays, the long non-Gaussian tails in 
the single particle freeze-out distributions are greatly reduced. In the longitudinal direction, 
where the non-Gaussian nature of the production point distribution is mostly defined by the 
system expansion dynamics (very similar for pions and kaons), the difference between the 
HBT radii and the R.M.S. values of freeze-out distribution is quite significant. The lambda 
parameter in the kaon case is generally larger than that for the pion case but due to a limited 
statistics exhibits large fiuctuations and is omitted from the figure. 

IV. CONCLUSIONS 

We have shown that the HBT radius parameters obtained from a Gaussian fit to the 
correlation function are equivalent to the the curvature of the relative freeze-out distributions 
at Ar = 0. Therefore, the fitted HBT radii are directly related to fitted Gaussian width of 
the two-particle position difference distribution near the origin. The lambda parameter is 
determined by the ratio of the value of the real relative distance distribution at the origin 
and the distribution chosen in a Gaussian form with the radii parameters determined by the 
fit to the shape of the correlation function. Studies with a dynamical model support the 
above statements. 
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FIG. 1. The single particle freeze-out position distributions and two-particle position difference 
distributions for pions with < pt < 200MeV/c and \y\ < 0.5 for r^j^e (perpendicular to the 
transverse momentum vector), rout (parallel to the transverse momentum vector), and ri^ng (along 
the beam axis). The Gaussian fit function is also shown in each case, along with the fit parameters. 
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FIG. 2. The fitted HBT radius parameters (solid circles), single-particle R.M.S. deviations 
(open squares), and the sigma of the Gaussian fit two-particle position difference distribution 
(open circles) as a function of pt for Rgide, Rout, Riong- The rapidity range is \y\ < 0.5. The bottom 
plot show the lambda parameter determined from the fit to the correlation function (solid circles) 
and from the two particle position difference distribution. 
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FIG. 3. The same as Fig. |2|, but for kaons. 



